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Abstract 

The magnetic properties of a quasi two-dimensional superconductor, Sn-metal graphite 
(MG), are studied using DC and AC magnetic susceptibility. Sn-MG has a unique 
layered structure where Sn metal layer is sandwiched between adjacent graphene 
sheets. This compound undergoes a superconducting transition at T c = 3.75 K at 
H = 0. The H-T diagram of Sn-MG is similar to that of a quasi two-dimensional 
superconductors. The phase boundaries of vortex liquid, vortex glass, and vortex 
lattice phase merge into a multicritical point located at T* = 3.4 K and H* = 40 
Oe. There are two irreversibility lines denoted by H g i (de Almeida-Thouless type) 
and Hgii (Gabay- Toulouse type), intersecting at Tq= 2.5 K and H' = 160 Oe. The 
nature of slow dynamic and nonlinearity of the vortex glass phase is studied. 

Key words: vortex glass phase, phase diagram, magnetic irreversibility 
PACS: 74.25.Ha, 74.25.Dw, 74.60.Ec, 74.80.Dm 



1 Introduction 

The existence of vortex glass phase in the mixed state of quasi two-dimensional 
(2D) high-T c superconductors such as YBa 2 Cu 3 7 in the presence of magnetic 
field (H) attracts considerable attention [1,2,3]. Miiller et al. [4] have pointed 
out for the first time the existence of the glassy phase in L^CuO^^Ba. They 
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have shown that the boundary between the vortex liquid phase and the vortex 
glass phase corresponds to the de Almeida-Thouless (AT) [5] line for the spin 
glass behavior. Since then there have been a number of experimental works 
on the vortex glass phase in quasi 2D superconductors. Experimentally it has 
been confirmed that the H-T diagram consists of the vortex glass, vortex 
lattice (Abrikosov lattice), and vortex liquid phases. The boundary between 
the vortex lattice phase and the vortex liquid phase (the line H a i) is of the 
first order, while the boundary between the vortex lattice phase and the vortex 
glass phase (the line H ag ) and the boundary between the vortex glass phase 
and the vortex liquid phase (the line H gl ) are of the second order. These 
boundaries merge into a multicritical point at T =T* and H = H*. 



Sn-metal graphite (MG) constitutes a novel class of materials having unique 
layered structures. This system can be prepared by reduction of an acceptor- 
type SnCl2 graphite intercalation compound (GIC) as a precursor material. 
Ideally, the staging structure of Sn-MG would be the same as that of SnCl2 
GIC. Sn layer is sandwiched between adjacent graphene sheets. These sand- 
wiched structures are periodically stacked along the c axis perpendicular to 
the basal plane of graphene sheets. Sn-MG undergoes a superconducting tran- 
sition at a critical temperature T c (= 3.75 K), which is close to that of bulk Sn 
[6]. The superconductivity occurs mainly in Sn metal layers, forming a quasi 
2D superconductor (type-II). 



In this paper we report results on the magnetic properties of Sn-MG from 
measurements of DC and AC magnetic susceptibility. Three methods are used 
to determine the H-T diagram: (i) the T dependence of the magnetization (M) 
in various states such as the FC (field-cooled) state, ZFC (zero-field cooled) 
state, IR (isothermal remnant ) state, and TR (thermoremnant) state, (ii) the 
T, H, h, and / dependence of the dispersion (Q[/ h or x') an d absorption (©'// h 
or x") °f the AC magnetic susceptibility, where h and / are the amplitude and 
frequency of the AC magnetic field, respectively, and (iii) the H dependence 
of M. We find that there are four lines (H ag , H g i, H g y and H a i) in the H- 
T phase diagram, separating vortex liquid, vortex glass, and vortex lattice 
phase. These lines merge at the multicritical point located at T* = 3.4 K and 
H* = 40 Oe. The lines H g i and H gV are irreversibility lines, intersecting at a 
critical point at T ' = 2.5 K and H' = 160 Oe. The line H gt has an AT-power 
law form, [5] while the line H g i> has a Gabay-Thouless (GT)-power law form 
[7]. The nature of possible slow dynamics and nonlinearity in the vortex glass 
phase is examined from the measurement of Q[/h and Q"/h as a function of 
u, h, T, and H. 
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2 Experimental procedure 



The method of sample preparation in Sn-MG is similar to that in other MG's 
previously reported [8,9,10,11]. SnCl 2 graphite intercalation compound (GIC) 
samples as a precursor material, were prepared by heating a mixture of natural 
graphite flakes (grade: RFL 99.9 S) from Kropfmiihl, Germany and an excess 
amount of SnCl 2 . The reaction was continued for three days at 400 °C in an 
ampoule filled with chlorine gas. The synthesis of Sn-MG was made by the 
reduction of SnCl 2 GIC. SnCl 2 GIC samples were kept for two days in a so- 
lution of lithium diphenylide in tetrahydrofuran (THF) at room temperature. 
Then the samples were filtered, rinsed by THF, and dried in air. The struc- 
ture of Sn-MG thus obtained was examined by (00L) x-ray diffraction, and 
bright field images and selected-area electron diffraction (SAED) (Hitachi H- 
800 transmission electron microscope) operated at 200 kV. No pristine SnCl 2 
and its hydrolisation products were left in the Sn layers. The details of sample 
preparation and characterization in Sn-MG will be reported elsewhere. 

The Sn-MG sample consists of many small flakes. Each flake has a well-defined 
c axis. The sample in the present work is regarded as a powdered sample 
with the c axis randomly distributed in all directions, since these flakes are 
randomly piled inside the sample capsule. The measurements of DC and AC 
magnetic susceptibility were carried out using a SQUID magnetometer (Quan- 
tum Design MPMS XL-5). Before setting up a sample at 298 K, a remnant 
magnetic field was reduced to less than 3 mOe using an ultra-low field ca- 
pability option. For convenience, hereafter this remnant field is noted as the 
state H = 0. (i)ZFC and FC magnetic susceptibility: The sample was cooled 
from 298 to 1.9 K at H — 0. After H was applied at 1.9 K, the zero-field 
cooled magnetization (M Z fc) was measured with increasing T from 1.9 to 6 
K. The sample was kept at 20 K for 20 minutes in the presence of H. Then 
the field cooled magnetization (M FC ) was measured with decreasing T from 
6 to 1.9 K. (ii) M-H loop: The sample was cooled from 298 K to T at H = 0. 
The measurement was carried out with varying H from to 3 kOe at T, from 
H = 3 to -3 kOe, and from H = — 3 to 3 kOe. (iii) M Z fc vs H: The sam- 
ple was cooled from 298 K to T at H = 0. The ZFC magnetization M ZFC 
at T was measured with increasing H from to 250 Oe. (iv) AC magnetic 
susceptibility. The dispersion (Q[/h) and the absorption (Q'(/h) was deter- 
mined from the in-phase and out-of phase first harmonic (uj) component in 
the AC magnetization measurement [12]. The dispersion (Q^/h) and the ab- 
sorption (Q'l/h) were described by <d[/h — x' + ^x^/A + 5x' 5 h 4 /8 + . . . and 
&[/h = X " + 3X3^74 + 5x^78 + • • • , where x ' 2n+1 and X 'L+i (n = l,2 r - •) 
are the real and imaginary parts of high-order nonlinear magnetic susceptibil- 
ity, respectively. In the limit of h = 0, Q'i/h and Q'[/h coincide with the linear 
susceptibility x' an d respectively. Before the measurement, the sample was 
cooled from 298 to 1.9 K at H = 0. (i) The T dependence of G'Jh and &{/h 
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Fig. 1. (a) - (f) T dependence of xfc and xzfc at low H (1 < H < 200 Oe) for 
Sn-MG. 

was measured with increasing T from 1.9 to 6 K in the presence of H . After 
the measurement, the sample was cooled from 6 to 1.9 K in the presence of 
the same H. Then the measurement was repeated at different H. (ii) The h 
dependence was carried out at fixed T by changing h (0.01 < h < 4 Oe) at 
/ = 1 Hz. (iii) The / dependence was carried out at fixed T in the presence 
of H by changing / (0.1 < / < 1000 Hz). 



3 Result 



3. 1 xfc and xzfc 



Figure 1 shows the T dependence of xfc an d Xzfc at l° w H. The susceptibility 
Xzfc deviates from xfc below a characteristic temperature for 1 < H < 200 
Oe, suggesting the irreversibility line in the H-T plane which separates the 
reversible region from the irreversible region. Figure 2 shows the T dependence 
of xfc and Xzfc ^or 1 < H < 3 kOe. For H > 300 Oe there is no noticeable 
difference between xfc and Xzfc at any T. The susceptibility (xzfc — Xfc) 
exhibits a very broad peak, which shifts to the low-T side as H increases. In 
the inset of Fig. 2(a) we plot these peak temperatures as a function of H. The 
value of H (= H p ) for the peak is much larger than that of a upper critical 
field H c2 (whose definition will be given in Sec. 4.1) at the same T, and it 
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Fig. 2. T dependence of \fc and xzfc at high H (1 < H < 3 kOe) for Sn-MG. The 
inset shows the plot of the peak temperature of xfc vs T as a function of H. The 
dotted line is a least-squares fitting curve. See the text for detail, (b) T dependence 
of xfc and xzfc at H = 3 kOe, which is a blow-up of Fig. 2(a). 



tends to zero around T c . This result suggests that an antiferromagnetic (AF) 
short-range order survives at high H where the superconductivity disappears. 
The origin may be due to the possible antiferromagnetism of nanographites 
in Sn-MG (see Sec. 4.3). 

Figure 3 shows the T dependence of the difference 5 (= Xfc — Xzfc) at various 
H. The difference 5 for H < H* (H* fx 40 Oe) decreases with increasing T 
and reduces to zero at a characteristic temperature. In contrast, the difference 
S for H* < H < 150 Oe decreases in two steps with increasing T, showing a 
local minimum and a local maximum. The difference S shows a small peak at 
2.25 K for H = 200 Oe and is almost equal to zero for any T for H = 250 Oe. 
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Fig. 3. T dependence of the difference 5 (= Xfc ~ Xzfc) at l° w H for Sn-MG. The 
solid lines are guides to the eyes. 

3.2 M IR and M TR 

In a mixed state of conventional type-II superconductors that have well-defined 
lower critical field H cl and upper critical field H c2 (H cl < H < H c2 ), the 
thermodynamically most stable state involves filling the system with flux- 
oids, exhibiting only a partial Meissner effect with the magnetic induction 
B = H + 4irM > 0. Even a small degree of inhomogenity is sufficient to 
pin some ffirxoids at local free energy minima, so that some fluxoids remain 
trapped in place where H is removed. In Sn-MG the values of M and B are 
strongly dependent on the magnetic states such as the ZFC, FC, IR (isother- 
mal remnant), and TR (thermoremnant) states. Experimentally, the values 
of M and B for these states are obtained as follows. First the sample was 
cooled from 298 to 1.9 K at H = 0. Then H was applied. The measurements 
of Mzfc an d Mj R were done with increasing T from 1.9 to 4 K. At each T, 
M Z fc was measured at H and then M IR was measured 100 sec later after 
the field was changed from H to Oe. Second, the sample was annealed at 
20 K for 1200 sec at H. The measurements of Mpc an d Mtr were done with 
decreasing T from 4 to 1.9 K. At each T, Mpc was measured at H =100 Oe 
and then M TR was measured 100 sec later after the field was changed from H 
to Oe. The values of M and B for each state thus obtained are described by 
Bzfc = H+A7rM ZFC , B ir = AttM ir , B fc = H+4irM FC , and B TR = AttM tr . 
Here B F q and Bzfc are the magnetic inductions in the FC and ZFC states 
with H, while Bj R and B TR are the magnetic inductions in the IR and TR 
states with H = 0, respectively. These magnetic inductions result from the flux 
trapping at pinning centers. Figure 4(a) shows the T dependence of M Z fc an d 
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Fig. 4. (a) T dependence of M Z fc, M ir , M fc , and M TR for Sn-MG. The definition 
of these magnetizations and the method of the measurement are described in the 
text. H = 100 Oe. (b) T dependence of 5 (= M FC - M Z fc), A (= M TR - M IR ), 
Mj R , and Mt R - (c) T dependence of (5 — A). The solid line is guide to the eyes. 

Mpc at H = 100 Oe (> if*), and M IR and M TR at H = 0. The T dependence 
of M/jj and Mtr is different from f?/# and B^ R by only a factor of 4-7T. The T 
dependence of M IR and M TR is similar to that of typical spin glass systems 
such as 3D Ising spin glass Fe c Mni_ c Ti0 3 [13], where the competition between 
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the nearest neighbor antiferromagnetic and ferromagnetic interactions leads 
to the spin frustration effect. The magnetization Mj R exhibits a broad peak 
around 3.0 K, while M FR decreases with increasing T. Both Mj R and Mtr 
reduce to zero at T c , implying the appearance of both B IR and B TR below T c . 
The deviation of M IR from M T r occurs below 3.1 K, showing the irreversible 
effect of magnetization. Note that the point at T = 2.94 K and H = 100 
Oe is located on the line H g i in the H-T diagram (see Sec. 4.1). Figure 4(b) 
shows the T dependence of A [=M TR — M IR = (B TR — B ir )/Ati} and 5 [= 
M FC — M Z fc — {Bfc — B zfc )/4:-k], as well as M IR and M TR . Figure 4(c) 
shows the T dependence of (5 — A). The difference (5 — A) is equal to zero 
at T above 3.25 K and exhibits a peak at 3.0 K where Mi R has a peak. It 
becomes negative below 2.7 K. The temperature region between 3.0 and 2.4 - 
2.7 K corresponds to the vortex glass phase for H = 100 Oe. 

Figures 5(a) - (c) show the T dependence of M FC) M ZFC) M ir , M TR) 5, A, 
and 5 — A for H = 30 Oe. The magnetization Mj R appears below T c , exhibiting 
a peak at 3.5 K. Note that a point at T = 3.5 K and H = 30 Oe is located 
on the line H a i in the H-T diagram (see Sec. 4.1). The deviation of Mj R from 
M T r occurs at T below 3.6 K, while the deviation of M Z fc from M FC occurs 
below 3.55 K. The difference (5 — A) exhibits a very sharp peak at 3.5 K, 
being equal to zero above 3.6 K and below 3.4 K. 



3.3 M vs H 



Figure 6 shows the hysteresis loop of the magnetization M at T — 1.9 and 5 
K. The M-H curve at 1.9 K shows a typical superconducting behavior with 
a negative local minimum around H = 160 Oe, corresponding to the phase 
boundary (the line H ag ) between the vortex lattice and vortex glass phase 
(see Sec. 4.1). The M-H curve is not reversible, partly because some fluxoids 
remain trapped at pinning centers when H is removed. The M-H curve at 
5 K is not perfectly proportional to H. The susceptibility Mj H is estimated 
to be 3.45 x 10~ 6 emu/g at H = 1 kOe, which is in contrast to the negative 
diamagnetic susceptibility of pristine graphite: x — ~ 30 x 10~ 6 emu/g for H 
along the c axis and —5.5 x 10~ 7 emu/g along the basal plane (so-called c 
plane) perpendicular to the c axis below 100 K [14]. 

Figures 7(a) and (b) show typical data of M Z fc vs H at various T. The sample 
was cooled from 298 K to T (< 4 K) at H = 0. The magnetization M ZFC at T 
was measured with increasing H (0 < H < 250 Oe). The magnetization M ZFC 
exhibits a single local minimum at a characteristic field for T < T c . Similar 
behavior has been reported in YBa2Cus06.92 by Shibata et al. [15]. It seems 
that a discontinuous change in M ZF c is observed just above the characteristic 
field for T* < T < T c , where T* = 3.4 K and T c = 3.75 K. The field of the 
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Fig. 5. (a) T dependence of M Z fc, M ir , M fc , and M TR for Sn-MG. H = 30 Oe. 
(b) T dependence of 5 (= M FC - M Z fc), A (= M TR - M IR ), M IR , and M TR . (c) 
T dependence of (5 — A). The solid line is guide to the eyes. 



local-minimum for M ZFC vs H, which is plotted as a function of T, forms the 
line H a i of the first order transition in the H-T plane (see Sec. 4.1). The jump 
height AMzfC) which is defined as the difference between the value of Mzfc 
at which Mzfc starts to be proportional to H and the local-minimum value of 
Mzfc vs H, rapidly decreases with increasing T and reduces to zero at T c . The 
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Fig. 6. Magnetization loop (M vs H) for Sn-MG. (a) T = 1.9 K and (b) 5 K. The 
method of the measurement is described in the text. The solid lines are guides to 
the eyes. 

jump height AM Z fc m ay be described by AM ZFC — M t — M a [16], where M t 
and M a denote the magnetization in the vortex liquid and the vortex lattice 
phase, respectively. The positive value of AMzfc implies that the vortex 
density in the vortex liquid phase is higher than that in the vortex lattice 
phase. The magnetization M ZF c at T = 3.7 and 3.8 K is almost proportional 
to H with a positive slope. The derivative dM ZFC /dH at 3.7 and 3.8 K shows 
a broad peak around H = 100 Oe, suggesting the existence of AF short-range 
order. The origin of the AF short-range order is due to nanographites in Sn- 
MG (see Sec. 4.3). Note that a discontinuous jump of M Z pc i n the vicinity 
of M ZF c = is an artifact due to the SQUID measurement. Such a jump in 
the magnetization is frequently observed when the magnetization changes its 
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Fig. 7. (a) and (b) Mzfc vs H at various T. The solid lines are guides to the eyes. 
The inset of Fig. 7(b) shows the jump in magnetization AAIzfc (which is defined 
in the text) SIS 3i function of T. 

sign. 

The lower critical field H c i(T) is defined as a field at which the curve of M 
vs H starts to deviate from the linear portion due to the penetration of the 
magnetic flux into the system: H c i(T = 1.9 K) = 20 ± 5 Oe. 

34 X' and X " 

Figure 8 shows the T dependence of the dispersion x' at various H, respec- 
tively, where / = 1 Hz and h — 0.1 Oe. The AC response to such a small 
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Fig. 8. T dependence of x' at various H for Sn-MG. h = 0.1 Oe. / = 1 Hz. The 
solid lines are guides to the eyes. 



Sn-MG 




Fig. 9. T dependence of x" at various H for Sn-MG. h = 0.1 Oe. / = 1 Hz. The 
solid lines are guides to the eyes. 

h is found to be linear in the whole range of temperatures and fields. The 
derivative dx'/dT shows a sharp peak at T c at H = 0, shifting to the low-T 
side with increasing H. The peak temperature for dx'/dT vs T plotted as a 
function of H in the H-T diagram forms the line H g i for H > H* and the line 
H al for H < H* (see Sec. 4.1). 

Figure 9 shows the T dependence of the absorption x" a t various H, where 
/ = 1 Hz and h — 0.1 Oe. The absorption x" at H = seems to have two 
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peaks at T c and 3.5 K. The peak at T c at H = shifts to the low-T side with 
increasing H . In contrast, the peak at 3.5 K at H = remains unshifted, 
becoming a shoulder even at 10 Oe, and tends to disappear above H = 30 
Oe. Note that the peak temperature of the main peak for x" y s T is slightly 
higher than that for d%' /dT vs T at the same H. The peak temperature for 
x" vs T is plotted as a function of H in the H-T diagram, forming the line 
H a i for H < H* and a new line H g \i for H > H* which deviates from the line 
H g i for H* < H < 150 Oe (see Sec. 4.1). The main peak of x" vs T is sharp 
for H < H* but becomes much broad for H > H*. 



3.5 G'Jh and Q"/h 



The AC response (Q[/h and Q"/h) for sufficiently large amplitude of the AC 
magnetic field becomes nonlinear below T c . In fact, as it will be shown later, 
the T dependence of Q"/h is strongly dependent on h for h > 0.2 Oe. This is 
in contrast to the case of h = 0.1 Oe, where the linearity is valid: Q[/h = x' 
and Q"/h = x" ■ Here a large amplitude of the AC magnetic field (h — 2 
Oe) was used to examine the effect of nonlinearity on the vortex glass phase. 
Figure 10 shows the T dependence of Q[/h and Q"/h for Sn-MG at various / 
(0.1 < f < 1000 Hz) in the presence of H (= 100 Oe), where h = 2 Oe. The 
T dependence of Q[/h is almost independent of /, while the T dependence of 
Q'l/h is strongly dependent on /. The shape of Q"/h vs T does not change 
with / but the basement of Q"/h proportionally increases with increasing / 
for 0.1 < / < 1000 Hz at any T below T c . As shown in Fig. 10(b) there are 
two broad peaks at 2.90 and 3.15 K for Q"/h with h = 2 Oe. These peaks 
do not shift with / for 0.1 < f < 1000 Hz within the experimental error in 
T (typically AT e = 0.01 K). No relaxation time effect is observed, suggesting 
out 1, where r is the relaxation time of the vortex system and uj (=2nf) is 
the measured angular frequency. This is very different from the case of typical 
spin-glass systems where ujt 1 and the peak of &[/h shifts to the high T 
side with increasing / [17]. We show that the T and / dependence of Q[/h 
and Q'l/h are explained in terms of the screening behavior arising from skin 
size effects. The absorption peak occurs when the electromagnetic penetration 
depth 8 S becomes comparable to the sample size L (see Sec. 4.2). The peak 
at 3.15 K coincides with that of x" at /i = 0.1 Oe and / = 1 Hz (see Fig. 9). 
The point at H = 100 Oe and T = 3.15 K is located on the line ti g v- On the 
other hand, the peak at 2.9 K is newly observed for h = 2 Oe. The point at 
H = 100 Oe and T = 2.9 K is located on the line H g i. Note that the data of 
dx' /dT vs T shows a peak at 2.94 K and d<5/dT has a local minimum at 3.0 
K. In this sense, the line H g i is the irreversibility line. 

Figure 11 shows the T dependence of Q"/h for Sn-MG at various / (0.1 < 
/ < 1000 Hz) in the presence of H (= 70 Oe), where h = 2 Oe. There are 
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Fig. 10. T dependence of (a) the dispersion Q[/h and (b) the absorption @ f [/h at 
various / (0.1 < / < 1000 Hz) for Sn-MG. H = 100 Oe. Because of the nonlinearlity, 
Q'l/h at h = 2 Oe is not equal to x" at /i = 0.1 Oe (see Fig. 12(b) for the nonlinearity 
at H = 0). The solid lines are guides to the eyes. 

two peaks at 3.1 - 3.15 K and 3.3 - 3.4 K. It seems that these peaks do not 
shift with /. This result suggests no occurrence of the relaxation time effect 
for 0.1 < f < 1000 Hz. The peak at 3.3 - 3.4 K coincides with that at 3.34 K 
of x" at /i = 0.1 Oe and / = 1 Hz (the data at H = 70 Oe are not shown in 
Fig. 9). The point at H = 70 Oe and T = 3.3 - 3.4 K is located on the line 
H g i<. On the other hand, the peak at 3.1 - 3.15 K is newly observed for h — 2 
Oe. The point at H = 70 Oe and T = 3.1 - 3.15 K is located on the line H g i. 
Figure 12 shows the h dependence of (a) Q[/h and (b) &{/h for Sn-MG at 
various T in the absence of H, where 0.01 < h < 4 Oe. Both Q[/h and Q"/h 
are not independent of h for h > 0.2 Oe, indicating the nonlinearity of the AC 
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Fig. 11. T dependence of &[/h at various / (0.1 < f < 1000 Hz) for Sn-MG. H = 75 
Oe. h = 2 Oe. The solid lines are guides to the eyes. 

response. The dispersion Q[/h linearly changes with h below T c , but is almost 
independent of h at 3.8 K just above T c . 



4 Discussion 

4-1 H-T phase diagram 

Our results obtained above are summarized in the H-T diagram. In Figs. 13(a) 
and (b), we make a plot of the peak temperatures of dx'/dT vs T and x" vs 
T, and the local-minimum temperatures of d<5/dT vs T as a function of H, 
and the local-minimum fields of Mzfc vs if as a function of T. We find that 
these data lie on the phase boundaries denoted by the lines H ag , H gh H g y, and 
H a i, which separate the vortex lattice, vortex glass, and vortex liquid phase. 
These lines merge at the multicritical point at T* = 3.4 K and H* = 40 Oe. 
The lines H ag , H g i, and H g y are of the second order, while the line H a i is of 
the first order. The line H C 2 is a crossover line separating the normal phase 
at higher temperatures from the vortex liquid phase, (i) On the line H g i, the 
data of dx'/dT vs T have a peak at fixed H, the data of dS/dT vs T have 
a local minimum, and the data of x" vs T have a peak for H > 160 Oe. (ii) 
On the line H g i', the data of x" vs T have a peak at fixed H . (iii) On the line 
H ag , the data of d<5/dT vs T have a local minimum at fixed H and the data of 
Mzfc vs H have a local minimum at fixed T. (iv) On the line H a i, the data 
of dx'/dT vs T and x" vs T have a peak at fixed H, the data of d<5/dT vs T 
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Fig. 12. h dependence of (a) ®Jh and (b) &[/h at various T for Sn-MG. / = 1 Hz. 
# = 0. 



have a local minimum at fixed H, and the data of Mzfc vs ^ have a local 
minimum at fixed T. (v) On the line H c2 , the data of Xzfc vs T come to take 
a saturated value Xzfc a ^ fixed H. The susceptibility Xzfc * s positive and 
weakly dependent on H. It takes 1.8 x 1CT 7 emu/g at H = 10 Oe, increasing 
with increasing H, and reaches a peak (= 4.63 x 10~ 6 emu/g) at H = 200 Oe. 
In turn it decreases with further increasing H, and reaches 3.75 x 10~ 6 emu/g 
at H = 1 kOe (see Fig. 3(a)). 

Using an empirical relation given by Werhammer et al. [18], H c2 (T = K) 
= — 0.69T c (dH c2 /dT)T c , the value of H c2 at T = K can be estimated as 
H&fi) = 1850 ± 20 Oe, where we use T c = 3.75 K and a slope {dH c2 / dT)rp c = 
—714.4 ± 6.8 (Oe/K) obtained from the linear relation in H c2 vs T in the 
vicinity of T = T c and H c2 = 0. The coherence length £(0) can be estimated 
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Sn-MG 




Fig. 13. (a) and (b) H-T diagram for Sn-MG. The peak temperatures of dx' /dT 
vs T (O), x" vs T (•), and the local-minimum temperatures of d<5/dT vs T (A) 
are plotted as a function of H. The local-minimum fields of Mzfc vs H (A) are 
plotted as a function of T. The vortex lattice (Abrikosov lattice), vortex liquid, and 
vortex glass phases are separated by the lines H a i, H ag , H g i, and H g v- These lines 
merge at the multicritical point located at (T* « 3.4 K, H* « 40 Oe). The line _ff C 2 
is a crossover line denoted by closed squares (■). The lower critical line H c \ is not 
shown in this figure. H C 2(T = OK) ~ 1.85 kOe. The solid lines for the lines H g i, 
H g [i , H ag , and H a i denote the least-squares fitting curves of the data for Hi vs T to 
Eq.(l) with i = gl,gl',ag, and al, respectively. The corresponding parameters are 
given in the text. 

as 420±20l using a relation H c2 (0) = $ /(2vr^(0) 2 ), where $ (= 2.0678 x 1(T 7 
Gauss cm 2 ) is a fluxoid. Note that our value of H c2 (0) for Sn-MG is on the 
same order as that of Sn films grown on an InSb (110) surface [19]. The value 
of H c2 (0) strongly depends on the thickness t of the Sn films: H c2 (0) = 3160 
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Oe for * = 420l and 2540 Oe for t = 2540A for H parallel to the Sn film. 



Here we assume that the T dependence of the line Hi (i = gl, gl', ag, and al) 
may be described by a power law form given by 

H = H*(l-T/T*) a(i) , (1) 



where a(i) is an exponent, and H* and T* are characteristic field and temper- 
ature, respectively. The least squares fit of the data of H vs T (40 < H < 250 
Oe) obtained from the peak temperature of dx'/dT vs T at H, to Eq.(l) with 
i = gl yields the parameters a(gl) = 1.43 ± 0.05, T* x = 3.90 ± 0.08 K, and 
Hgi = 730±20 Oe. The exponent a(gl) is close to that (= 1.50) predicted from 
the AT theory [5], suggesting that the line H g i is the so-called irreversibility line 
separating a higher-temperature vortex liquid phase from a lower-temperature 
vortex glass phase. The least squares fit of the data of H vs T (70 < H < 160 
Oe) obtained from the peak temperature of x" vs T at H, to Eq.(l) with 
i = gl' yields the parameters a(gl') = 0.57 ± 0.10, T* v = 3.57 ± 0.09 K, and 
H* v = 330 ± 10 Oe. The exponent a(gl') is close to that (= 0.50) predicted 
from the Gabay-Thoules (GT) theory [7]. Such an occurrence of the AT-GT 
crossover has been reported in granular high T c superconductors [20,21]. No 
GT-power law form has been observed in a single crystal of YBa 2 Cu 3 7 with- 
out granularity [20] . These results suggest that the AT-GT crossover behavior 
is indicative of the frustrated nature of the system due to disorder and granu- 
larity. This crossover is an analogy to the Ising-Heisenberg crossover behavior 
observed in Heisenberg-like spin glasses [22]. Such a crossover occurs when 
the magnetic field collapses the random local anisotropy field. The AT and 
GT lines in SG systems represent the longitudinal and transverse freezing, 
respectively. 

Similarly, the least-squares fit of the data of H vs T (40 < H < 150 Oe) 
obtained from the peak temperature of dx'/dT vs T at H, to Eq.(l) with 
i = ag yields the parameters a(ag) = 2.06 ± 0.05, T* g = 4.77 ± 0.80 K, and 
H* g = 435 ±50 Oe. Note that the lower critical field H c i(0) is much lower than 
H* . The least-squares fit of the data of H vs T (0 < H < 40 Oe) obtained 
from the peak temperature of dx'/dT vs T at H, to Eq.(l) with i = al yields 
the parameters a(al) = 1.02 ± 0.03, T* t = 3.70 ± 0.02 K, and H* al = 580 ± 50 
Oe. The value of a(al) is in good agreement with the mean-field value (= 1) 
derived from the Ginzburg-Landau theory [23], where the correlation length 
varies with (1 - T/T^)" 1 / 2 . 

The H-T diagram of Sn-MG is similar to that of a quasi-2D superconductors 
(type-II). The field for the multicritical point (= if*) is described by the 
dimensional crossover field that separates the region of 2D thermal fluctuations 
(for H > H*) and the 3D thermal fluctuations (H < H*). In this sense, the 
H-T diagram may be universal for any quasi-2D superconductors when the H 
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axis and T axis are scaled by H/H* and T/T c , respectively. For YBa2Cu307 
with T c = 90 K, the multicritical point is located at T* = 76 K and H* = 90 
kOe [2]. The ratio T*/T c for YBa 2 Cu 3 07 is 0.84, which is almost equal to that 
(= 0.91) for Sn-MG with T* = 3.4 K and T c = 3.75 K. 

4-2 Nature of the lines H gi and H gV 

Here we discuss the origin of the irreversibility lines H g i> and H g i for H* < H < 
160 Oe. We find the following results on the T dependence of AC magnetic 
susceptibility at small AC field {h — 0.1 Oe) and large AC field (h = 2 Oe). 
(i) The linear absorption x" & t h = 0.1 Oe exhibits a peak only on the line 
Hgi/. (ii) The nonlinear absorption Q'[/h exhibits a peak on the line H g i< and 
a peak on the line H g \. 

It is assumed that there is an irreversibility line T g i{H) (or the line H g i(T)) 
in the H-T plane which separates an Ohmic region [T > T g i(H)] from an 
non-Ohmic region [T < T g i(H)}. Above T g i(H), a voltage V is proportional to 
current /, whereas below T g i(H), V shows an extremely nonlinear dependence 
where the linear resistivity is zero and the electric field E is described by 
E m exp(— A/ J* 1 ) [24]. Here A and p are positive constants and J is the 
current density. 

Why does the peak of the absorption for small AC field amplitude appear on 
only the line H g y for H > H*7 Why does the peak of the absorption for large 
AC field amplitude appear on both the lines H g i and H g ii for H > H*l Similar 
problem was addressed first by Geshkenbein et al. [24] and later by Steel and 
Graybeal [25] . According to their skin size effect hypothesis, the peak of ©'// h 
appears when the electromagnetic penetration depth 5 S = (c 2 p/ f) 1 ^ 2 /2tt is on 
the order of the system size L, or when u coincides with the characteristic 
angular frequency uj peak given by 

iu peak = 0.8c 2 p(T,H)/L 2 , (2) 

where p(T,H) is the electrical resistivity. Then Q"/h is a function of only 
oo/uj peak , and has a peak at uj — ou peak . When 5 S becomes infinity, the field 
penetrates the system completely and 0'///i tends to zero. In the opposite 
limiting case (S s — > 0), the screening is complete and Q'{/h = 0. 

First we consider the peak temperature T peak (H) of Q"/h at any h, where 
Tpeak(H) > T g i(H). We assume that the response of the system is linear 
(Ohmic). The current density is of the order of J pea k = ch/(4nL). The resis- 
tivity p(T,H) = E(J peak )/ J P eak, which is independent of J pea k, decreases with 
decreasing T and tends to zero at T gl (H). The peak temperature T peak (H) 



19 



can be determined from the matching condition (uj = u pea k) given by Eq.(2): 
T = T^ ak for A = 0.1 Hz and at T = T^ ] ak for f 2 = 1 kHz. Then we have 

a relation p{T^ ak , H) / p(Tpl ak , H) = 10 4 . Since there is no shift of T pea k(H) in 
Q'l/h with / (0.1 < / < 1000 Hz) in the present work, it follows that the dif- 
ference of peak temperatures defined by AT = T^ ak — T^ ak should be less than 
0.01 K. Because of the sharp drop in resistivity with decreasing T, T peak (H) 
thus obtained is very close to but above T gl (H). The value of T peak (H) is in- 
dependent of h because Eq. (2) is independent of h. The resistivity is finite at 
T = T peak (H): p min = p(T peak ,H) = 1.25uL 2 /c 2 . In summary, Q'(/h vs T has 
a peak at T peak (H) just above T g i(H), which is independent of h. 

Next we consider the peak temperature T peak (H) of Q"/h at any h, where 
Tpeak(H) < Tgi(H). We assume that the response of the system is highly non- 
linear (non-Ohmic). The effective resistivity is defined by p e ff(J peak ,T, H) = 
E(J pe ak)/ Jpeak, where J P eak = ch/(4nL). As a function of T, p e ff decreases 
with decreasing T and tends to zero at a temperature below T g i(H). As a 
function of J pea k, Peff has a local minimum at J pea k = (Ap) 1 ^, and is equal 
to zero at both J pea k = and J pea k = oo. When J pea k is smaller than (Ap) 1 ^ 
[or h < 4:7rL(ApY^ /c], p e ff (Jpeak) decreases with decreasing J pea k (or with 
decreasing h). The peak temperature T peak (H) can be estimated from Eq.(2) 
with p(T,H) = p eff (Jp eak ): p e ff(J P eak) = 1.25cuL 2 /c 2 at T = T peak (H). The 
value of J P eak (or h) is uniquely determined because of the relation between 
Jpeak and h. The peak temperature T pea k(H), which is lower than T g i(H), 
increases with decreasing h and approaches the line T g i(H) from the low-T 
side. In summary, Q'(/h vs T has a peak at T peak (H) just below T g i(H) for 
sufficiently large h satisfying the above condition. 

The lines H g i and H g y intersect at a critical point located at Tq = 2.5 K and 
Hq = 160 Oe. The line H g i line for H > H' Q no longer separates the irreversible 
region of the vortex glass phase from the reversible region of the vortex liquid 
phase. In this sense, the vortex glass phase may connect smoothly to the vortex 
liquid phase at the critical point at Tq and H' Q [3]. As shown in Fig. 3, in fact 
the T dependence of 5 at H > 200 Oe is rather different from that at 150 Oe. 
The difference 5 at H > 250 Oe is almost equal to zero at any T, indicating 
the reversibility of magnetization on the line H g i(T). 

The line H g i(T) may correspond to the AT line of spin glasses [5]. Then it is 
expected that the relaxation time of Sn-MG, r, diverges as a power law form 
r = T [T/T g i(H) — l]~ x [1], as T approaches the irreversibility line T g i(H) from 
the high-T side. Here r is a characteristic relaxation time and a; is a dynamic 
critical exponent. Experimentally, no relaxation time effect is observed in the 
T and / dependence of Q[/h and Q'(/h on the line H g i(T) (see Sec. 3.5). This 
result indicates that the condition ujt 1 is satisfied, in spite of crossing the 
line H g i(T). Thus the power law form is not the case of the present system. 
What is the origin of r satisfying the condition lut <^ 1? The absorption of the 
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AC field is mainly governed by the hopping of the vortices over the barriers 
between different metastable states. The depinning time r for such intervalley 
transition is obtained from a characteristic relaxation time r by an Arrhenius 
form t r exp(U /kBT), where U is the activation barrier energy between 
the different metastable states [24,25]. Then r smoothly changes as T changes 
and wr < 1 on the line H g i(T). 

In summary, it is concluded from the experimental results described above 
that the properties of the vortex glass phase in Sn-MG are characterized as 
follows, although most of them are not sufficiently understood at present, (i) 
The magnetization M Z fc deviates from M FC below T gi (H) for H > H*. (ii) 
The magnetization Mir deviates from Mtr below T g i(H). The T dependence 
of Mir and Mtr is very similar to that of Mzfc and Mpc in typical spin 
glass systems, respectively [17]. (iii) The difference (5 — A) takes a positive 
large value for the vortex glass phase, while it is equal to zero for the vortex 
liquid phase and it takes a negative value in the vortex lattice phase, (iv) The 
relaxation time r shows no behavior of divergence on the line T g i(H). (v) The 
non-Ohmic behavior is observed in the vortex glass phase, while the Ohmic 
behavior is observed in the vortex liquid phase. 



4-3 Antiferromagnetic short-range order 

There are several results supporting the existence of the AF short-range order. 
The first evidence is the T dependence of the peak field H p as shown in the 
inset of Fig. 2. The field H p is much higher than H c2 at the same T: H p m 3.2 
kOe and H c2 ~ 490 Oe at 3.0 K. The field H p increases with decreasing T and 
may be described by a power law form given by H p = H*{1 — T/T*)*, where 
j3 is an exponent, and H* and T* are characteristic field and temperature, 
respectively. The least squares fit of the data of H p vs T (1 < H < 3 kOe) to 
the power law form yields the parameters f3 = 0.41 ±0.02, T* = 3.82 ±0.02 K, 
and H* = 5.95±0.05 kOe. The exponent (3 is close to /3 — 1/2 for the molecular 
field theory. The temperature T* is a little higher than T c . The second evidence 
is the magnitude of the DC magnetic susceptibility: x — 3-45 x 10~ 6 emu/g 
at H = 1 kOe and T = 5 K. This value of x is positive and is very different 
from that of pristine graphite which is diamagnetic (see Sec. 3.3) [14]. 

What is the origin of the AF short-range order? The existence of the AF 
short-range order has been observed in other MG's such as Bi-MG [26]. In 
these MG's the metal layers would generate internal stress inside the graphite 
lattice, leading to the break up of adjacent graphene sheets into nanographites. 
Fujita et al. [27] and Wakabayashi et al. [28] have theoretically suggested that 
the electronic structures of finite-size graphene sheets depend crucially on the 
shape of their edges (zigzag type and armchair type edges). Finite graphite 
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systems having zigzag edges exhibit a special edge state. The corresponding 
energy bands are almost flat at the Fermi energy, thereby giving a sharp peak 
in the density of states (DOS) at the Fermi energy. As a result, conduction 
electrons are localized around the zig-zag edges of nanographites, forming lo- 
calized magnetic moments. Recently Harigaya [29] has theoretically predicted 
that the magnetism in nanographites with zigzag edge sites depends on the 
stacking sequence of nanographites. The antiferromagntic spin alignment be- 
comes energetically favorable when the adjacent graphene sheets have a A-B 
type stacking sequence along the c axis, like pristine graphite [30]. The AF 
short-range order in nanographites may coexist with the superconductivity 
occurring in Sn sheets even at low T and low H . 



5 CONCLUSION 

Sn-MG is a quasi-2D superconductor, showing a superconducting transition 
at T c = 3.75 K for H = 0. The H-T diagram of Sn-MG consists of four 
lines (H ag , H g i, H g y, and H a i) separating vortex lattice phase, vortex glass 
phase, and vortex liquid phase. These lines merge into the multicritical point 
located at T* = 3.4 K and H* = 40 Oe. The irreversibility lines H gt and H gV 
intersect at a critical point at Tg = 2.5 K and H' Q = 160 Oe. The line H gi 
has an AT-power law form, while the line H g y has a GT-power law form. The 
slow dynamics and nonlinearity in the vortex glass phase is an analogy to 
that in the spin glass phase observed in 3D Ising spin glasses. Further studies 
will be required to understanding the nature of the vortex glass phase from 
measurements of resistivity (the nonlinear relation of current-voltage) and the 
absorption and dispersion at frequencies (/ ^> 1 kHz). 
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